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Abstract — Transparency  of  the  four  atmospheric  extremely  high  fre¬ 
quency  (EHF)  window  ranges  located  around  35,  90,  140,  and  220  GHz  is 
obscured  by  precipitation  (rain,  wet  snow),  by  suspended  particles  (fog. 
cloud,  haze,  dust),  and  by  water  vapor.  An  assessment  is  made  of  the 
quantitative  picture  (i.e.,  models  versus  experiments  and  theory),  upon 
which  estimations  for  general  radio  path  behavior  can  be  based.  Useful 
models  are  provided  for  calculating  attenuation  based  upon  measurable 
meteorological  variables:  rain  rate,  liquid  water  content,  humidity,  tem¬ 
perature,  and  pressure.  Information  currently  available  is  not  yet  com¬ 
plete  enough  to  provide  accurate  predictions  under  all  climatological 
conditions.  Emphasis  is  on  recent  advances  in  formulating  the  physical 
basis  for  modeling  transparency  and  on  a  discussion  of  some  of  the 
principal  remaining  uncertainties. 
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I.  INTRODUCTION 

New  requirements  for  high  data-rate  communication  links,  for 
active  and  passive  sensors,  and  in  radio  astronomy  have  kindled 
growing  activity  in  the  extremely  high  frequency  (EHF)  (30  to 
300  GHz)  region  of  the  radio  spectrum.  Progress  is  spurred  by 
available  bandwidth  and  size  reduction,  and  by  advantages  over 
infrared  and  optical  systems  in  penetrating  an  opaque  atmosphere. 
Nevertheless,  atmospheric  propagation  limitations  dominate  most 
considerations  in  the  advancement  of  applications-adverse 
weather  causing  attenuation  by  rain,  wet  snow,  and  fog  and  by 
water  vapor,  is  of  serious  consequence  to  EHF  system  perform¬ 
ance.  In  recent  years,  numerous  research  and  measurement  pro¬ 
grams  have  been  carried  out  to  make  these  detrimental  effects 
predictable.  A  wealth  of  material  has  been  accumulated  and  was 
carefully  reviewed  [l]-[9].  Here  we  limit  the  discussion  to  at¬ 
tenuation  in  transparent  portions  of  the  EHF  range,  somewhat 
loosely  defined  by  a  transmission  factor 


T  =  exp  (-0.2303 A)  >0.0 1  or  A  <20  dB.  (1) 

Expression  (1)  implies  that  at  least  one  percent  of  the  trans¬ 
mitted  power  can  be  received.  The  total  radio  path  attenuation  A 
in  decibels  is  a  measure  of  the  energy  extracted  from  a  plane  wave 
by  the  atmospheric  propagation  medium.  For  a  given  situation, 
attenuation  A  can  be  evaluated  when  the  number  ol^all  types 
of  absorbers  and  scattcrers  within  a  radio  path  is  known.  A 
specific  attenuation  a  is  introduced  when  the  mean  absorber 
density  is  constant  along  a  path  length  /.,  yielding 


TABLE  I 

EHF  WINDOW  FREQUENCY  RANGES  OE  THE  ATMOSPHERE 
(SEE  FIG,  1) 


WINDOW  FREQUENCY  RANGE  CLEAR  AIR  ATTENUATION 

(Sea  Level,  15°C,  50  RH) 


ID 

f 

• 

GHz 

.1  , 

dB/km 

Wl 

25 

to 

50 

0.1 

to  0.5 

M2 

70 

to 

115 

0.5 

to  2 

M3 

125 

to 

160 

2 

to  5 

MA 

200 

to 

250 

5 

to  1 0 

In  clear  air  one  can  expect  relatively  low  values  for  a  (i.e.. 
useful  lengths  L  for  a  given  value  of  A )  in  four  distinct  frequency 
ranges  (see  Table  I),  for  which  the  metaphor  “window"  is  used. 
The  frame  for  these  windows  is  set  up  by  molecular  resonance 
absorption  of  H20(centered  at  22.  183.  and  325  GHz)  and02 
(centered  at  60  and  119  GHz)  as  illustrated  in  Fig.  1  (10].  It 
appears  that  about  150  GHz  of  the  EHF  band  are  available  for 
signal  transmission.  A  closer  look  reveals  that  the  windows  are 
veiled  by  water  vapor  and  suspended  hydrometeors  and  darkened 
by  precipitation. 

Recently  Crane  [11],  [1]  compared  all  readily  available  clear 
air  EHF  window  data  with  model  calculations  for  a  standard  con¬ 
dition  (sea  level,  20°C,  7.5  g/m3  H20).  Attenuatiim  a,  within 
measurement  uncertainty,  was  supported  by  about  130  data 
points  from  26  references;  and  for  a  zenith  path  attenuation  A,. 
starting  at  standard  sea  level  conditions,  about  60  data  points 
from  22  sources  were  available.  He  agrees  with  previous  workers 
that  an  empirical  correction  is  required  to  reconcile  observations 
wilh  H20  line  calculations  and  that  information  is  missing  to 
establish  the  physical  basis  for  such  a  correction.  This  mini¬ 
review  goes  one  step  further  in  that  it  is  assumed  that  a  theoretical 
model  modified  by  an  empirical  term  [10]  correctly  renders  the 
window  attenuation  dependences  and  is  then  reduced  to  power 
law  expressions  in  terms  of  meteorological  variables. 

II.  TRANSMISSION  PRINCIPLES  AND  EHF  CLIMATOLOGY 

Analysis  of  T  (I)  is  not  trivial.  The  loss  of  plane  wave  power 
by  the  cumulative  effects  of  all  absorbers  in  a  radio  path  of 
length  /,  is  tallied  by  the  total  attenuation 
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a  =  <4//,  dB/km.  (2) 
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A  I  |aw(x)  3  a„.(.v)  F  o,  (v)  +  o/;(.v)|  r/.v  dB.  (3) 

■'(i 

wheic  a  is  the  specific  aitenuation  of  an  air  volume  clement  over 
a  path  increment  t/.v  and  the  subscripts  R.  W.  V.  D  designate  the 
absorbers  rain,  suspended  hydrometeors  (fog.  cloud,  haze,  water 
cluster),  water  vapor,  and  dry  air.  respectively. 

The  physical  state  of  air  at  a  point  v  along  the  path  is  identified 
by  five  basic  meteorological  parameters  pOv,  w.  and  R  described 
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Fig.  1.  Specific  attenuation  at  sea  level  over  the  frequency  range  1-250  GHz  for  various  relative  humidities  (0  to  100  percent), 

including  fog  (0.1  g/m3)  and  rain  (R  =  10  mm/h). 


and  ranged  below. 

A.  Moist  air: 

p(x)  =  (P-e)  kPa,  (4a) 

dry  air  pressure  where  Pis  the  barometric  pressure; 

0(x)  =  300/7,  (4b) 

relative  inverse  temperature  (1.5  to  0.95  for  T  =  200  to  3 1 5  K); 

u(x)  =  7.2I9e0  g/m3,  (4c) 

water  vapor  concentration  where  e  is  the  partial  vapor  pressure  - 
or  by 

RH(x)=41.51e/(0sX  10(10-9  8340))  «  29<?018 

<100  percent,  (4d) 

relative  humidity  (0=0,  is  the  dew  point  at  relative  humidity 
(RH)  =  100  percent). 

Water  vapor  variability  at  sea  level  (p  =  101  kPa,  0  =  1 .016  or 
22°C)  is  typically 


Dry 

Normal 

Humid 

Saturated 

1 

10 

17 

20 

g/m3 

RH  =  5 

50 

85 

100 

percent 

B.  Suspended  hydrometeors  are  described  by  the  liquid  water 
concentration  w,  which  relates  approximately  to  optical  (0.55 
pm)  visibility  (/(km)  [4|, 

w(x)  =  0.01 1/(7'  7  g/m3.  (4e) 

A  schematic  categorization  can  be  made  by 


Convective 


Cluster  (2 1 1 

Haze 

Fog  1 44 | 

Stratus 

Cloud 

tv  *  I0~3 

10-2 

10-1 

1 

5  g/m3 

lia  17 

l.l 

0.27 

0.07 

0.03  km 

C.  Precipitation  originates  as  a  highly  statistical  event  within 
clouds  suspended  in  saturated  air.  Its  vertical  distribution  is 
separated  into  two  regions  by  the  height  of  the  0°C  isotherm, 
which  can  vary  between  ^6  km  and  ground  level  depending  on 
season  and  latitude.  The  lower  part  is  mostly  liquid  drops,  and 
the  upper  region  consists  of  frozen  particles  with  occasional 
supercooled  droplet-loadings  by  strong  updrafts. 

Yearly  statistics  on  local  point  rain  rates  /?(*,-)  have  proven 
useful  in  modeling  rain-induced  attenuation  effects  [13]— {18] , 
[43],  [45).  Rain  rate  can  be  related  to  the  percent  time  iK,  a 
given  value  occurs  over  the  period  of  an  “average"  year:  the  effec¬ 
tive  rain  cell  extent  Lr/L\  and  the  instantaneous  suspended 
liquid  water  concentration, 

wR  =  mR(x)  g/m3.  (40 

In  terms  of  these  variables,  a  typical  local  rain  may  be  classified  as 
follows  (horizontal  path.iL  =  10  kni): 


Drizzle 

Steady 

Heavy 

Downpour 

Cloudburst 

R  =  1 

5 

20 

100 

250 

mm/h 

<R  ~2 

0.5 

0.07 

0.001 

0.0001 

percent/yr 

LrH  -  1 

1 

0.7 

0.35 

0.2 

m  *0.1 

0.07 

0.05 

0.04 

0.04 

The  simple  coefficient  scheme  reveals  some  fundamentals  of 
rain  that  require  consideration  in  attenuation  models.  Changes  in 
the  factor  m  indicate  rain  rate-dependent  characteristics  of  drop 
size  distributions.  Widespread  steady  rain  occurs  more  uniformly 
and  favors  small  drop  sizes  (<l  mm  diameter)  which  stay  in  the 
air  longer.1  Heavy  showers  are  more  localized,  favor  larger  drops 
and  occur  less  frequently. 

1  Fall  velocities  in  still  air  range  between  1  and  9  m/s  for  drop  dia¬ 
meters  between  0.2  and  7  mm. 
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Fig.  2.  Normalized  (1  g/m3, 300  K)  specific  H20  attenuation  over  the  frequency  range  5-1000  GHz  for  five  phase  states:  water 
1),  rain  2),  suspended  hydrometeors  3),  moist  air  4)  (numbers  shown  indicate  selected  H20  attenuation  peaks),  and  pure 
water  vapor  5)  [10] . 


Three  snow  classes  may  be  approximated: 


Dry 

Moist 

Wet 

R$  =  0.1 

1 

5 

mm/h 

T  =  -10 

-5 

0 

♦5<,C, 

where  a  snow-equivalent  rate  Rs  can  be  obtained  from  a  heated 
rain  gauge. 

Climatological  extremes  for  locations  in  the  continental  U.S. 
span  the  following  range. 


extremes  of  continuum-type  water  (curve  1)  and  resonance-type 
vapor  (curve  5)  spectra.  Parameterization  of  specific  attenuation 
on  the  basis  of  (4)  is  treated  in  the  following  manner. 

•  Rain  attenuation  can  be  fitted  to  (13) 

aR  =  oRx  dB/km,  (5) 

bypassing  elaborate,  lengthy  Mie  calculations  which  require  drop 
shape,  size  distributions,  and  the  complex  dielectric  constant  of 
water  [12] -[14].  Drop  diameters  and  EHF  wavelengths  (0.1-5 
mm)  are  roughly  comparable,  thus  causing  appreciable  attenuation 


Rain  (R) 

Water  Vapor  (u) 

Temperature  (fl) 

Annual  Rate 

Greatest  Monthly  Concentration 

Mean  Annual  Relative  Temperature 

Minimum 

180  mm,  Phoenix 

6  5  g/m3,  Denver 

1.042  (14.8°C),  San  Francisco 

Maximum 

1620  mm,  New  Orleans 

21  g/m3,  Miami 

0.984  (31.7°C),  Phoenix 

At  this  point  it  is  assumed  that  for  a  given  radio  path  the  geo-  via  resonance  (Mie)  scattering.  The  frequency-dependent  coef- 

metric,  climatological,  and  temporal  statistics  of  the  meteoro-  ficient  a  and  exponent  x  (Table  II)  were  calculated  using  drop- 

logical  measurables  (4)  are  available  and  that  models  and/or  s‘ze  spectra  of  Laws  and  Parsons. 

theories  exist  to  correlate  (pOv  +  w  +  R)  combinations  with  A  Considerable  effort  has  been  made  to  measure  and  to  describe 
(3).  rain  attenuation.  Rain  cell  extent  LR  and  path-averaged  rain  rate 

R  are  obtained  from  cumulative  distributions  of  attenuation  A  or 
III.  SPECIFIC  ATTENUATION  point  rate  R  at  many  locations  over  long  time  periods,  preferably 

In  the  EHF  window  ranges,  it  is  possible  to  formulate  simpli-  several  years.  Statistical  estimation  methods  and  models  for  the 
fied  expressions  for  the  specific  attenuation  o.  Water  in  both  quantity  (aL)R  are  discussed  in  detail  by  Ippolito  [2],  [3],  A 

liquid  and  vapor  states  obscures  window  transparencies  for  large  bibliography  exists  on  the  topic  of  computing  aR  distri- 

ground-based  system  applications.  Frequency  dependence  and  butions  for  general  radio  paths  on  the  basis  of  point  rainfall 

magnitude  of  H20 attenuation  is  distinctly  different  for  the  same  observations.  More  recent  models  are  those  by  Lin  or  Dutton 

absorber  thickness  1  mm  (later  see  (13))  in  water,  rain,  suspended  et  a!.  (U.  S.)  [15],  [43] ;  Morita  (Asia)  [17] ;  Misme-Waldteufel 

droplet,  moist  air,  and  pure  water  vapor  forms,  as  depicted  in  Fig.  (Europe)  [42] ;  and  Crane  (worldwide)  [18] .  A  typical  require- 

2.  The  three  atmospheric  phase  states  lie  somewhat  between  the  ment  in  system  design  is  to  know  for  a  given  location  the  par- 
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TABLE  II 

COEFFICIENTS  a-d  AND  EXPONENTS*, y  FOR  CALCULATING  SPECIFIC  ATTENUATION  or  (3)  IN  THE  FOUR 

EHF  window  Ranges  (see  fig.  l) 


EHF 

Window 

Frequency 

f 

Rain  (Low) 

R  <50 

a  x 

Rain  (High) 

R  =  25-200 

a  x 

Snow 

*S  *S 

Haze.  Fog,  Cloud 

b  y 

Water  Vapor 
c  c* 

Dry  Air 

d 

GHz 

Reference 

[131 

(191,1201 

110) 

X  10  —  7 

X  10—7 

25 

0.11 

1.08 

0.14 

1.02  ( 

0.05  dry  1.7 

0.30 

7.0 

0.017 

0.44 

0.022 

Wl 

35 

0.24 

1.02 

0.34 

0.91  < 

0.25  moist  .9 

0.59 

6.4 

0.010 

0.27 

0.035 

50 

0.48 

0.91 

0.66 

0.82  ' 

0.5  wet  .6 

1.17 

5.9 

0.017 

0.44 

0.270 

70 

0.80 

0.81 

0.86 

0.79  / 

0.3  dry  1.5 

2.20 

4.8 

0.032 

0.83 

0.336 

W2 

90 

1.03 

0.76 

0.94 

0.78  ^ 

1.4  moist  1 

3.45 

3.6 

0.049 

1.27 

0.054 

115 

1.12 

0.73 

0.98 

0.77  l 

2.5  wet  7 

5.17 

2.2 

0.082 

2.11 

0.230 

130|,110| 

120 

1.16 

0.72 

0.99 

0.77  i 

0.5  dry  1.3 

5.53 

1.9 

0.087 

2.22 

0.800 

W3 

140 

1.25 

0.71 

1.00 

0.77  < 

2.5  moist  1 

6.95 

0.9 

0.116 

2.98 

0.029 

160 

1.30 

0.70 

1.01 

0.76  ' 

4.0  wet  .9 

8.36 

-0.1 

0.190 

4.88 

0.024 

200 

1.48 

0.66 

1.06 

0.76  ( 

[  0.7  dry  1 

10.9 

-1.6 

0.28 

7.14 

0.024 

W4 

220 

1.47 

0.66 

1.05 

0.76  \ 

3.0  moist  1 

12.1 

-2.1 

0.24 

6.20 

0.026 

250 

1.47 

0.66 

1.03 

0.76  ' 

[5.2  wet  1 

13.4 

-2.8 

0.30 

7.63 

0.029 

AVERAGE  OF  REPORTED  FIELD  DATA 

a  Reference 

(moist)  »s  (we0  Reference 

b 

Reference 

c 

Reference 

Wl 

35 

0.2 

0.4(1)  1 7  ] 

0.4(1) 

[7| 

0.009(3) 

1 7 ) , [11] 

W2 

90 

0.9(1)* 

[7| 

1.8(5) 

3.0(6)  J7),  1 191 

4.5(5) 

(71 

0.04(2)  [4],  [7],  |11] 

W3 

140 

1.1(3) 

171 

3.0(8) 

5(1) 

m,  [i9i 

6(1) 

[71 

0.13(5)  [4], (7), |11) 

W4 

220 

15(8) 

[201 

3.5(10) 

6(1) 

|20) 

10(2) 

1 20  J 

0.3(1)  |4|,|201,[11) 

*  Digits  in  parentheses  give  the  standard  devaition  from  the  mean  in  terms  of  the  final  listed  digits, 
a  (db/km)/(mm/h). 
b,  c  (dB/km)/(g/m3). 
c'  (dB/km)/percent. 
d  (dB/km)/kPa. 

ticular  rain  rate  (attenuation)  that  is  not  exceeded  by  more  confirmed,  would  add  very  little  to  aw  since  even  relatively  high 
than  0.01  percent  (1  h)  of  the  time  over  a  year’s  period.  Costly  predicted  number  densities  (<  1013  cm-3)  [21]  only  yield  low 
design  decisions  hinge  on  this  rain  information.  For  example,  concentrations,  w  <  I0~ 3  g/m3. 


the  1  h/yr  rain  rate  expectation  in  the  continental  U.  S.  has 
extremes  between  7  (Arizona  desert)  and  115  mm/h  (south 
Florida).  Rain  attenuation  statistics  for  a  terrestrial  path  in  mid- 
U.  S.  are  shown  in  Fig.  3,  as  calculated  with  Crane's  prediction 
model  (zone  D2)  [18],  extended  in  frequency  beyond  100  GHz 
by  results  from  [  1 3 J . 

•  Snow  is  difficult  to  model  since  it  is  highly  inhomogeneous. 
Around  0°C,  flakes  carry  varying  amounts  of  free  water.  During  a 
snow  storm  n  the  northeast  of  the  U.  S.,  the  equivalent  liquid 
concentration  was  found  to  vary  between  0.2  and  0.8  g/cm3  [  1 9) . 
Data  on  and  R$  are  few  and  have  been  fitted  to  (5).  The  re¬ 
sulting  parameters  «sand.xs(seeTable  II)  yield  only  first  estimates. 

•  Suspended  hydrometeor  attenuation  is  formulated  by 


Interrelations  between  the  “EHF  weather”  parameters  (4)  are 
not  treated,  although  H20  is  an  incessantly  phase-changing  agent 
in  the  hydrological  cycle  between  condensation  sites  and  the 
atmospheric  water  vapor  content.  For  example,  at  high  relative 
humidities  (RH  >  90  percent),  a  small  portion  of  water  vapor  is 
converted  into  submicron  droplets  by  water  uptake  of  aerosol 
particles  leading  to  the  formation  of  haze  (w  <  10“ 2  g/m3)  [22] . 
Heavy  fog  and  clouds  can  exceed  0.1  g/m3,  causing  appreciable 
specific  attenuation  a^.  Cloud  attenuation  in  W3  and  W4  can 
rival  rain  attenuation  depending  upon  the  extent  Lw  (0.1  to  3  km) 
of  the  absorption  cell. 

•  Water  vapor  attenuation  is  modeled  by 

aF  =  c(p/\0\ )vOy  =  c'(p/ IQ\ )  RH/01 7  dB/km  (7) 


aw  =  bw$y  dB/km. 


^  where  y  =  1 .8,  2.4,  and  3.2  in  Wl ,  W2,  and  W3,  4,  respectively. 


The  calculation  of  a*  is  based  upon  the  Rayleigh  absorption 
approximation  of  Mie  scattering  [14],  Diameters  of  suspended 
particles  are  below  0.1  mm  and  any  detailed  drop-size  distribu¬ 
tion  reduces  to  a  liquid  water  concentration  w.  The  Debye  model 
for  the  complex  dielectric  constant  of  water  enters  into  the 
computation  of  the  parameters  A  and)’  (see  Table  II).  Unlike  rain. 


The  attenuation  depends  on  absolute  (u)  humidity  limited  by  RI1 
(note  the  steep  temperature  dependence  (4d)),  and  varying 
greatly  with  both  time  and  location  [23] .  On  a  hot  (35°C)  sum¬ 
mer  day,  the  vapor  concentration  can  reach  V  =  40  g/m3  at  RH  > 
90  percent. 

The  scale  factor  c  has  its  origin  in  the  absorption  linesol  ll20 


a  strong  temperature  dependence  can  be  noticed.  Neutral  water  vapor.  In  principle,  it  can  be  evaluated  from  a  line  by-line  sum- 
dusters  (five  to  30  molecules),  whose  existence  is  not  directly  mat  ion  of  all  line  features,  each  contributing  at  a  given  frequency 


IEEE  TRANSACTIONS  ON  ANTENNAS  AND  PROPAGATION,  VOL.  AP  31.  NO.  1,  JANUARY  1983 


131 


ATTENUATION  -  dB 

Fig.  3.  EHF  window  rain  attenuation  predictions  for  a  terrestrial  ( h  = 
0  km)  path  of  1  km  length  using  for  a  mid-U.S.A.  region  the  Crane 
model  [  18 1 . 


TABLE  III 

WATER  VAPOR  CONTlNUUM:COEFFIClENTS  DUE  TO  LINES  [27 1  ABOVE  1  THZ 


Equation 

k 

X 

y 

z 

Line  Shape 

(8a) 

0.0014 

1.0 

2 

2.51 

MGR 

(GR)  x  f/vQ 

(Modified  GR) 

OS 

(8b) 

0.0064 

1.05 

1.4 

2.05 

VW 

Van  Vleck-Weisskopf 

(24), [10] 

w 

(8c) 

0.012 

1.05 

1.4 

2.04 

FL  =  GR 

full  Lorentz  or  Gross 

[24] 

£ 

(8d) 

0.064 

1.05 

1.2 

1.17 

SL 

single  Lorentz 

(10) 

(8e) 

0.11 

1.1 

1.2 

1.10 

MVW 

(VW)  x  V/ 

(Modified  VW) 

z 

w 

a 

(80 

0.049 

1 

2.1 

2.0 

VW  I 

for  the  local 

[8|,  [11|, Table II 

g 

(8g) 

0.061 

1 

2.1 

1.22 

GR  » 

lines 

[29[, [30| 

3 

as  prescribed  by  a  line  shape  function  [24] .  Five H20 line  bases,  LB1  to  5,  have  been  reported: 


Line  Base  ID 

LB1 

LB2 

LB3 

LB4 

LB5 

Number  of  lines 

38350 

17201 

2277 

151 

30  ♦  ac 

Highest  Frequency,  THz 

537 

126 

15 

3 

i 

Remark 

AFGL  tape 

revised 

rotational 

radio 

EHF  approximation 

AFGL 

spectrum 

astronomy 

Reference 

|25] 

1261 

1 27 1 

1 28 ) 

1 10| 

The  low  frequency  wings  of  about  180  stronger  (e.g„  6  X  I04 
dB/km  at  v0  =  2.64  THz,  p  =  101  kPa,  v  =  1  g/m3,  300  K)  air- 
broadened  lines  below  20  THz  make  a  contribution  in  the  EHF 
range  [4] .  For  LBS,  30  local  lines  are  formulated  explicitly,  and 
far-wing  contributions  from  lines  centered  above  1  THz  are  fitted 
below  400  GHz  to  a  continuum  (slowly  varying  /-dependence) 
absorption 

ac=k{pl\Q\)if6y(fl\0OY  dB/km.  (8) 

A  theoretical  fit  to  (8)  was  accomplished  on  the  basis  of  LB3  by 
subtracting  the  local  lines  below  I  THz  from  all  2277.  Applying 


to  each  line  five  different  standard  shape  functions  led  to  the 
results  summarized  in  Table  ill.  Since  the  scale  factor  k  varies 
over  a  wide  range  (e  g..  1:70  at  100  GHz),  the  obvious  question 
is  “which  shape  function  is  correct?"  One  criterion  is  the  fact 
that  integrated  line  absorption  has  to  be  finite.  A  shape  function 
/■'(/).  normalized  to  unity  at  the  line  center  frequency  p0.  should 
have  a  constant  area  ny{y  is  line  width)  under  the  curve  when  / 
varies  between  0  and  °°.  A  test  can  be  defined  by 

fi  =  Id/TT-y)  /  F(f)df\  1=0  (0) 

Jo 
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and  used  for  cheeking  the  various  line  shapes.  The  most  simple 
function  is  a  single  Lorentzian  (SL)  which  in  normalized  form 
reads  as 

^SL  =a2bj[a2  +(1  —b)2]  (10) 

where  a  =  yli>0  and  b  -  //u0 .  Integration  of  (9)  in  closed  form  is 
not  possible  for  all  shape  functions.  A  numerical  integration  over 
a  limited  range  b  =  10“  3  to  103  and  a  =  0.001  (u0  =  3000,  7  = 

3  GHz)  with  five  choices  of  F(J )  (see  Table  III  and  |36] )  results 
in  various  amounts  of  deviation  from  (9): 

S  =  0.0000  (FL  =  GR),  0.0019  (SL).  0.0038  (MVW). 

0.0082  (MGR),  0.64  (VW). 

Although  departures  from  (9)  are  caused  predominantly  by  incor¬ 
rect  high  frequency  far-wing  contributions,  they  also  make  the  re¬ 
spective  low  frequency  far-wing  contributions  suspect.  Without 
getting  involved  in  discussing  merits  of  the  various  shape  function, 
we  note  that  only  FL  and  GR  oblige  rigorously  to  (9).  The  MGR 
and  VW  shapes,  based  on  (9)  as  well  as  other  arguments,  are  not  ap¬ 
plicable  to  molecular  line  absorption  over  a  wide  spectrum  range 
(0.01  to  20  THz)  and  their  use  for  the  water  vapor  spectrum 
should  be  discontinued  [32] .  More  refined  line  shape  treatments, 
possibly  capable  of  predicting  seven  orders  of  magnitude  in  at¬ 
tenuation  over  many  octaves  require  additional  parameters  [32 ]  - 
[34 1  besides  i>0  and  7. 

In  fitting  experimental  water  vapor  EHF  absorption,  many 
workers  applied  the  MGR -shape-with  the  result  that  large  deficits 
were  left  unaccounted  for  and  then  labeled  "excess"  [4],  [35  J . 
The  question  whether  the  excess  is  a  line  shape  problem  or  is 
caused  by  unaccounted  absorbers  (e.g.,  dimers  and  clusters  of 
H20)  has  been  a  long-standing  source  for  speculations  in  the 
correct  modeling  of  water  vapor  attenuation  [4] ,  [7] ,  [8[ ,  [11], 
[29],  [35].  Equation  (8)  is,  in  the  window  ranges,  the  major 
(up  to  85  percent  at  100  GHz)  contributor  to  av  (7).  The 
empirical  continuum  (80,  which  was  used  in  [8] -[1 1 1 ,  provides 
too  high  an  estimate  for  the  scale  factors  c  (Table  II)  and  cz 
(Table  V)  in  ranges  W3  and  W4,  leaving  (8g)  the  better  choice. 
This  Rice-Ade  continuum  (see  Table  III)  is  supported  by  theory 
assuming  the  Lorentzian  shape  (10),  that  is  (8g)  at  (8d). 

•  Dry  air  attenuation  mainly  due  to  oxygen,  can  be  estimated 
with  the  relation 

aD  =  dpO3  dB/km.  (11) 

Oxygen  makes  a  small,  though  always  present,  contribution.  The 
coefficient  d  (Table  11)  was  calculated  with  the  model  detailed  in 
[10],  The  line-by-line  summation  includes  overlap  effects  of  the 
60  GHz  absorption  band  and  nonresonant  02  plus  pressure- 
induced  Nj  absorption. 

•  Smoke  and  suspended  dust  do  not  pose  serious  EHF  attenua¬ 
tion  problems  [20].  A  visibility  in  smoke  of  U  as  0.2  km  or  a 
suspended  mass  concentration  (dust)  of  0.01  g/m3  exhibit 
specific  attenuations  on  the  order  of  0.15,  0.5,  1 ,  and  3  dB/km  in 
Wl  to  W4,  respectively. 

•  In  summary,  the  dominant  role  ofHjO  in  influencing  atmos¬ 
pheric  transparency  (1)  for  terrestrial  applications  is  evident. 
Concentrations  of  the  various  phase  states  are  highly  variable. 
Rain  presents  the  most  serious  limitation  to  system  perform¬ 
ance.  Assessing  individual  contributions  to  specific  attenua¬ 
tion  in  (3)  is  difficult  from  measurements  of  A.  Reliability,  pre¬ 
cision,  and  scale  of  supporting  data  (4)  compromise  the  quality  of 
most  observations.  Progress  lies  more  in  controlled  laboratory 
experiments,  where  it  is  possible  to  study;  elements  (except  rain) 
of  the  accounting  sum  (3)  in  isolation. 


Specific  attenuations  due  toll-,0, (5)  to(K).  increase  with  fre¬ 
quency:  also,  temperature  plays  a  role.  At  220  Gllz.  the  effective¬ 
ness  to  absorb  has  increased  with  reference  to  35  Gil/  by  different 
factors  17=  [a(220)/a(35)] 

Absorber/  Rain  Snow  Suspended  Hydrometeor  Moist  air 

n  5  10  16  24 

IV.  RADIO  PATH  ATTENUATION 

Attenuation  (3)  is  evaluated  from  [pOv  +  w  +  R)  distributions 
along  a  radio  path.  For  a  short  horizontal  path  a  set  of  average 
conditions  (4)  can  be  assumed  and  specific  attenuations  are 
simply  multiplied  by  the  geometric  path  length  /. . 

The  cumulative  attenuation  of  a  zenith  radio  path  requires 
height  profiles  of  p  -  0  ~  v,  u\  and  R  as  prescribed  by  in  sim 
data  or  by  a  synthetic  atmosphere  (two-dimensional  model). 
Liquid  water  content  W,  total  precipitable  water  vapor  F,  rain 
cell  depth  /,#  (15|-[18|.  and  dry  air  attenuation  dz  determine 
zenith  attenuation  which  is  modeled  by  [36] 

Az  =aR/.R  +bzW  +  ctV  +  dz  dB.  (12) 

where  the  water  and  vapor  column  densities  are 

W=  I  w(h)dlt  and  V-  I  Vji)dh  mm.  (13) 
'0  '0 

Both  quantities  If  and  V  can  be  measured  simultaneously  with  a 
ground-based  microwave  radiometer  over  time  scales  from 
minutes  to  months  [38] ,  [39]  .The  coefficients  b.,  c.,  and  dz  have 
been  evaluated  for  arctic,  standard,  and  subtropical  atmospheres 
specified  in  Table  IV.  Excess  radio  path  length  l.t:  due  to  refrac¬ 
tive  delay  provides  additional  information  about  the  air  mass; 
dispersive  effects  from  molecular  resonances  [10]  are  neglected. 
Results  of  three  sets  of  air-mass-specific  coefficients  are  listed  in 
Table  V.  The  example  of  A.  in  Fig.  4  is  for  the  U.  S.  Standard 
Atmosphere  [37].  Calculation  of  A.  was  performed  “layer-by- 
layer”  in  a  spherically  stratified  atmosphere  encompassing  48 
slabs  between  h  =  0  and  30  km.  The  ratio  b./cz  demonstrates  the 
contribution  of  suspended  hydrometeors  to  window  absorption, 
relative  to  water  vapor. 

Slanted  radio  paths  in  clear  air  with  an  elevation  angle  0  >  6° 
against  horizontal  follow  a  secant  law  [40]  .  i.e.. 

Ai>-Az/sm<t>  dB.  (14) 

For  modeling  through  rain,  clouds,  and  low  angle  situations  see 
[2],  [15]-[17] ,  [45]. 

V  CONCLUSION 

Atmospheric  EHF  attenuation  is  a  function  of  the  radio  path 
absorber  population  subject  to  considerable  variations.  Model¬ 
ing  runs  the  gamut  from  meteorological  observables,  through 
specific  attenuations,  to  theoretical  and  experimental  aspects 
of  window  transparencies.  Results  front  a  wide  spectrum  of 
recent  work  have  been  reviewed.  Intricate  “tnolccule-by -mole¬ 
cule.”  “line-by-line,"  "drop-by-drop"  and  “layer-by-layer"  cal¬ 
culation  procedures  have  been  reduced  to  simple  expressions 
allowing  rapid  estimations  of  window  transparencies.  No  at¬ 
tempt  was  made  to  cover  all  possible  angles  of  the  subject 
matter.  Rather,  a  direct  line  of  thought  was  followed  in  order 

to  convev  some  imderstandino  <>l‘  ITIt'  ■•11,,,..  . 
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TABLE  IV 

SELECTED  PROPERTIES  OF  THREE  MODEL  ATMOSPHERES 


ATMOSPHERE 


( RH=1 00%> 


AVERAGE 

REFRACTIVITY 

EXCESS 

TEMPERATURE 

path  length 

U.S.Std 

Subtrop 


10 7  N(h)dh 


TABLE  V 

COEFFICIENTS  bz,  Cz,  AND  dz  FOR  CALCULATING  ZENITH  ATTENUATION  A z  IN  EHF  WINDOW  RANGES 
(SEE  FIG.  4)  FOR  MODEL  ATMOSPHERES  SPECIFIED  IN  TABLE  IV 


EHF  Center 

Window  Frequency 

Model 

Prediction 

Experiments 

Atmosphere 

ID 

bz 

[10! 

«z 

d* 

b*/c  * 

5z 

dz 

Reference 

GHz 

dB/mm 

dB 

dB/mm 

dB 

X10-2 

xlO-2 

Arctic  1 

1.71 

1.10 

0.24 

155 

W1 

35 

US.  Standard  2 

1.40 

0.91 

0.20 

154 

0.5(3) 

0.16(2) 

|4 1 . 1 1 1 1 

Subtropic  3 

1.00 

0.84 

0.19 

119 

1 

4.49 

5.90 

0.36 

76 

W2 

90 

2 

4.83 

4.67 

0.32 

103 

2.6(5) 

0.3005) 

|4|.|U| 

3 

4.63 

4.19 

0.30 

111 

[30|,|10I 

1 

S.41 

14.6 

0.19 

37 

W3 

140 

2 

6.58 

10.8 

0.17 

61 

9(2) 

0.2(1) 

HI .  1 1 1 1 .  |4 1 1 

3 

7.53 

10.1 

0.16 

75 

1 

5.92 

32.0 

0.16 

19 

W4 

220 

2 

7.75 

25.7 

0.15 

30 

15(5) 

0.10) 

HMIII.HU 

3 

10.3 

23.3 

0.14 

44 

2401) 

1 30 1 
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U.  S.  Std.  Atm.  76 


Fig.  4.  One-way  zenith  attenuation /t2  through  the  U.S.  Standard  Atmos¬ 
phere  [ 37 ]  over  the  frequency  range  5-300  GHz  in  75  GHz  steps  for 
dry  1),  moist  2),  humid  3),  and  saturated  air  4)  and  for  4)  containing 
a  rain-bearing  cloud  5). 


Curve 

ID 

Delay 

LE 

Vapor 

V 

Suspended  Water 

W 

m 

mm 

mm 

1 

2.283 

0 

0 

2 

2.302 

2.87 

0 

3 

2.378 

14.36 

0 

4 

2.473 

28.71 

0 

5 

28.71 

1.18 

tion  effects.  The  outcome  of  observations  can  be  predicted, 
provided  the  physical  variables  (4)  are  available.  In  general,  agree¬ 
ment  with  experimental  data  is  within  limits  of  their  stated  un¬ 
certainty  (typically  better  than  ±  15  percent).  The  theoretical 
basis  for  water  vapor  continuum  absorption  was  discussed  in  more 
detail  but  still  awaits  confirmation  of  a  correct  and  practical  line 
shape  function.  Presently,  the  number  and  accuracy  of  experi¬ 
mental  data  decreases  sharply  above  100  GHz.  There  is  a  future  in 
gathering  more  accurate  data  suitable  to  make  adjustments  to  the 
proposed  coefficient  scheme. 
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